Resonance energy and charge pumping through quantum SINIS contacts 
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We propose a mechanism of quantum pumping mediated by the spectral flow in a voltage-biased 
SINIS quantum junction and realized via the sequential closing of the minigaps in the energy spec- 
trum in resonance with the Josephson frequency. We show that the dc current exhibits giant peaks 
at rational voltages. 

PACS numbers: 73.23.-b, 74.45. +c, 74.78.Na 
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Quantum pumping in mcsoscopic structures offers 
a unique possibility for studying and direct manipu- 
lating the fundamental quantum characteristics of the 
nanoscale objects (see, for example, 0, 0, Q and refer- 
ences therein). Most of the previously proposed quan- 
tum pumps (including the superconductor-quantum-dot 
junction of Ref. |^) used the cyclic adiabatic processes in 
order to avoid the relaxation which was believed to smear 
out the quantum behavior In this Letter we pro- 
pose and investigate a realization of the nonequilibrium 
resonance charge pump in a form of a superconductor- 
normal-superconductor (SNS) junction which essentially 
relies on both the discrete level dynamics due to the adi- 
abatic variation of scattering parameters of the SN con- 
tacts and on the relaxation processes in the continuum 
quasiparticle spectrum. We discuss an exemplary de- 
vice, the symmetric voltage-biased SINIS structure (see 
Fig. Pi consisting of two superconducting leads (S) cou- 
pled via the (tunnel) barriers (I) and the quantum ballis- 
tic normal conductor (N). The dc current (the pumped 
charge) exhibits giant peaks at rational bias voltages pro- 
vided the chemical potential of the normal conductor 
is varied in certain compliance with the Josephson fre- 
quency. This requires a nonequilibrium electron distri- 
bution achieved under the condition that the Josephson 
frequency is higher than the inelastic relaxation rate but 
smaller than the Andreev levels spacing in the normal 
conductor. The relaxation of this distribution is of crit- 
ical importance since it enhances the pumping efficiency 
to the degree that the dc current peaks greatly exceed 
the rectification current (Shapiro steps) observed in the 
equilibrium state of the same junction, thus providing un- 
ambiguous manifestation of the quantum pumping effect 
(see the discussion in Ref. Q). 

The normal conductor of the length d > v x /\A\ is cho- 
sen to have a single conducting mode; here v x is the Fermi 
velocity of the conducting mode and |A| is the supercon- 
ducting gap in the leads. The energy spectrum of sub-gap 
Andreev states, e n {<j>), has a large number N ~ \A\d/hv x 
of levels each being a function of the phase difference (f> 
between the superconducting leads; a typical spectrum 
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FIG. 1: Charge pump through a quantum channel with a 
variable chemical potential. The bias voltage is V = Vi, — Vr. 



is shown in Fig.JSfa). For a finite strength Z of barriers 
between the leads and the normal conductor, the levels 
are separated by the minigaps Ae2n at <f> = 2irk and try 
minigaps Ae w at <j> = 7r(l + 2k), where k is an integer |5|. 
Magnitudes of these minigaps depend on the interference 
phase a' = k x d + 8 between the waves incident on and 
reflected from the barriers (k x is the Fermi wave vector 
of the mode and S is the scattering phase). The mini- 
gaps Ae w disappear in the resonance determined by the 
condition sin a' = 0. In its turn, Ac2tt disappear in the 
anti-resonance where sin a' = ±1. The Fermi velocity 
and thus the interference phase a' can be tuned by the 
gate voltage V g as shown in Fig. Thus one can close 
either the gaps at <f> = 2wk or the gaps at <j> — ir(\ + 2k) 
sequentially adjusting a proper time dependence of V g . 
Choosing the bias voltage V — (h/2e)d<fi/dt such that 
<t> = at the time when the gap Ae27r is closed and (j) = tt 
at the time when Ae n is closed, and so on, a resonance 
energy and charge pumping via spectral flow of excita- 
tions from states below — |A| to states above +|A| and 
back will take place realizing an "Archimedean screw" Q 
in energy space. The particles moving upwards will reach 
continuum at e = +|A| with the distribution correspond- 
ing to the equilibrium at e = — |A| and vice versa. The 
subsequent relaxation of "wrong" distributions in contin- 
uum is accompanied by large dissipation leading to a dc 
current component under a dc bias voltage. It is essential 
that the gaps disappeared sequentially: vanishing of just 
one set of gaps either at 2?rk or at 7r(l + 2k) does not 
lead to the energy or charge pumping. A possible way to 



2 



i e/|A| (b) 




FIG. 2: Energy spectra Eq. © for a SINIS contact with 
Z = 0.5 and Ad/hv x — 10. (a) Dark lines: In the reso- 
nance sin a' = 0; light (red online) lines: in the anti-resonance 
| sin a' | = 1. (b) Pathes connecting states below — |A| to those 
above +|A| for trajectories m — n = of Fig. |3 



realize this device is to bridge two superconductors by a 
carbon nanotube subject to a gate voltage. A similar 
technology was used in Ref. Q to design semiconductor 
nanowires with tunable supercurrent. 

Resonances.- To close the mini-gaps in resonance with 
<f> one needs to tune the amplitude and frequency of the 
gate voltage in accordance with the bias voltage. Let us 
plot in Fig. |3 two sets of points in the <f), a' plane: (A) 
points 4> — 2irk, a' = 7r(fc + 1/2) (shown by crosses) and 
(B) points 4> — ir(2k + 1), a 1 — irk (circles). Resonance is 
achieved when, changing a' and <j> with time, one obtains 
a trajectory that passes first through a point A, then 
through a point B, then again through A, and so on. 
For linear trajectories parametrized as 2a' — uj g t + 2irk, 
(j) — 7r = ujjt, where ujj — 2eV/H, this requires ujj/uj g = 
(l+2m)/(l+2n). Some of these trajectories are shown in 
Fig. by dashed lines. The resonance trajectories open 
continuous paths in the e, <j) space connecting states below 
— |A| to those above |A|, see Fig. Hfb). The gate voltage 
can be varied within finite intervals along the trajectories 
[broken lines in Fig. [3] equivalent to the straight lines 
since sin 2 a' = sin 2 uj g t/2. More practical is to apply a 
sinusoidal gate voltage with a frequency O. For given m 
and n, the trajectory passing through the same points is 



is the Pauli matrix in Nambu space, and 



A 
A* 



The superconducting gap is A = |A|e ±Ic ^ 2 for x > d/2 
and x < —d/2, respectively, while A = for —d/2 < x < 
d/2. For simplicity we model the normal reflections at 
the interfaces as being produced by (5-function barriers 
U(x) = I5(x - d/2) + I5(x + d/2). 

In the normal region the particle, e ±lq+x , and hole, 
e^ iq ~ x , waves have amplitudes and , respectively. 
The upper or lower signs refer to the waves propagating 
to the right = (u + , v~) or to the left ■0 < = v + ). 
The particle (hole) momentum is q± = k x ±e/hv x . Scat- 
tering at the right and left barrier couples the amplitudes 
of incident and reflected waves la : 



1>r = s ^ , i>l = S L i>< ; S - 



S21 S22 



(3) 



The scattering matrices for the right and left barriers 
are S R = S(x2) and S L — S(\i), respectively, where 
Xi = -~4>/2 while \2 — 0/2. Components of the S matrix 
for 5- like barriers and energies |e| < |A| are |!j 



S\i{x, Z) = S 2 2(x, -z) = - 
S 12 (x,Z)e-^ = S2i( X ,zyx 



(U 2 -V 2 )(Z 2 + iZ) 
U 2 + (U 2 - V 2 )Z 2 
UV 

~ U 2 + {U 2 - V 2 )Z 2 



Here Z = ml/h 2 k x is the barrier strength and U = 
2- 1 / 2 [l + i^\A\ 2 - e 2 /e] 1/2 , V = U*. The waves at dif- 
ferent ends of the normal channel have different phases 
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where a — k x d, [3 — ed/hv x . Using unitarity of S ma- 
trix, S^S — 1, the condition of solvability of Eqs. © and 
Q assumes a compact form 



ISnpsin 2 a' 



+ |5 12 | 2 cos 2 (0/2) = sin 2 (/3 + 7 ) . (5) 



(ed/hv x )V g = [(1 + 2n)7r/2] sin(Ot) (1) e 2 ^ 



Here a' = a + 5 and the scattering phase 8 is intro- 
duced through cot S = Z; the phase 7 is defined as 



where </> — 7r = ujjt with ojj = 20(1 + 2m); the equivalent 
u) g = 20(1 + 2n). The amplitude of the gate voltage 
should be odd rational of irhv x /2ed; note that eV g <C A 
for long contacts with N ^> 1. 

Spectrum.— We use the Bogoliubov-deGennes equations 



2m dx 2 



E F + U(x) 



e"°> — Sn/S22- Equation (J5J determines the energy 
spectrum of a SINIS contact. Such spectrum has been 
extensively studied by many authors (see e.g. [l(ij]). 

In what follows we focus on long contacts, d ^> 
fit>a,/|A|, with a large number of levels N. Examples 
of the spectra are shown in Fig. E[a). All gaps at 
4> = 7r(l + 2fc) disappear at the same resonance parameter 
sin a' = 0. Similarly, all gaps at (j) — 2nk disappear when 
I sin a' I = 1. This follows from Eq. Q due to unitarity 
a z ip + HiJj = etp , (2) l^ii | 2 + l^il 2 = 1. The low-energy levels, e„ < A, have 



the form e„ = ±£q + 7rfrv x n/d where 



to explicitly demonstrate the resonance properties of the 
sub-gap energy spectrum of a SINIS structure. Here a z 



hv 7 



arcsm 



T 2 cos 2 — 
2 



(1 - T 2 ) sin 2 a' , (6) 
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FIG. 3: Linear resonance trajectories. 



n is an integer, T = (1 + 2Z 2 ) 1 , and we use 7< 1. The 
spectrum has energy gaps at <j) = n(l + 2k) 



Ae^ = (2Hv x /d) arcsin | sin a' | yl — T 2 



The gaps Ae2?r at = 27rfc are given by the same expres- 
sion with I sin a' I replaced with | cosa'|. For a transpar- 
ent contact, T = 1, minigaps disappear. 

Current- For time-dependent E g and (j) the current 
can be calculated by varying Eqs. J2J with respect 
to 8 A and the chemical potential of the normal conductor 
controlled by the gate voltage: 



6E F = SE g {x) = 



-eSV g , \x\ < d/2 
0, \x\>d/2 



Note that variations of j3, 7, and of the scattering phase 
5 resulting from changes in k x are negligible. If the mag- 
nitude |A| is constant, 8A = iS\A, we find 



(h/2e) (SxrJtiR ~ SxlJul) = 5e n - Q n SV & 



g > 



(7) 



where Q n = e J (\u n \ 



2 )dx is the charge on level 



n; the integral in Q n is extended only over the length of 
the normal conductor since \u n \ 2 — \v„\ 2 — in an open 
superconducting region. We denote 



-{ieh/2m) [u* n (du n /dx) + v^(dv n /dx) 



JnR,L are the values of J n near the right (left) barrier. 

Let e n , x> an d E g adiabatically depend on time as a 
parameter. Summing up the Andreev states having the 
distribution function /„ we find the energy balance 



g tne 
■ HI 



2/» 



(8) 



/sg _ _|_ jsg-j j2 j g current carried by the sub- 
gap states, (h/2)dxR,L/dt = — sVr^l = ±eV/2 due to 
symmetry; Vr^l and I^ L are the potentials and currents 
in the leads. In Eq. JHJ, the work done by the bias voltage 
VI SS and by the gate voltage (the first term in the r.h.s.) 
produce the energy change of the Andreev states. 



Keeping the phases fixed in Eq. Q we find Q n — 
de n /dV g . Since e„ is a function of two parameters, V g 
and 4>, Eq. ||SJ contains only the partial derivative 



' ft ^ dcj> { 



2/n) 



(9) 



where dcf>/dt = 2eV/H was used. This equation for mally 
coincides with the known time- independent result [illl2|. 

Note that the average currents satisfy I g = Ir—Il = 
where I g is the gate current (see Fig. ^| . Therefore I is 
the average current through each lead. It is determined 
by the longer time of the particle drift through the levels 
while the average gate current has also contributions due 
to much more rapid jumps of the distribution function 
each time when one of the levels joins the continuum, 
thus7 g ^-£„Q n (l-2/„). 

Charge pumping- Consider the simplest case where 
the gate voltage varies along the broken line in Fig. 
equivalent to a straight line m = n = 0. Let the deviation 
from the resonances be 



uijt , 2a = ujjt + 2aa 



(10) 



where ao is determined by initial conditions, and calcu- 
late the average current as a function of ao. As the phase 
point moves along the trajectory with ao ^ 0, the dis- 
tance between nearest levels does not become less than 
Ae. For low energies from Eq. JBJ) and small ao, 

Ae = (fti^/d) ^47-2(1 - T 2 )a 2 +ojj(t - t ) 2 . 

The probability of Zcncr tunnelling between the levels is 

p = exp [~2ttv x T 2 (1 - T 2 )al/ujjd] 

for all <j> s» irk. The Andreev states Eq. (JSJ are well- 
defined under an adiabatic condition that p is small far 
from the avoided crossings, u>j -C v x T 2 (l — T 2 )/d. 

Consider first an ideal adiabatic process when transi- 
tions between the levels are absent and the distribution 
/„ on each level is constant. Equation (JSJ) yields 

BR r)f 1 

sm[2((3 + l)]f € ^ = ^\Si2? sin(0 - n) . (11) 

Note that d/3/de = d/hv x . The derivatives dj/de ~ 
d\S\ 2 /de - (A 2 - e 2 )- 1 / 2 if Z ~ 1; they can be neglected 
for long contacts if e is not very close to A. We put <f> — 
(f)' + 7r in the average (de n /d(j)) = (27r) _1 J q ™ (de n /d(j)) d(j) 
and see that sin[2(/3 + 7)] in Eq. (fTT|l does not change 
its sign, while de/dcf> does if the particle stays on the 
same level. The average current vanishes for any ao =/= 0: 
Pumping is absent without interlevel transitions followed 
by relaxation of the electronic distribution in continuum. 

If p is finite, the particles move up or down the spec- 
trum while the distribution relaxes in continuum due to 
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fast escape through the potential barriers. The latter 
requires wj <C v x T/d which is not as strict as the adia- 
batic condition above. The average current is now finite 
for nonzero ag. When p reaches p = 1 for a = 0, a 
free drift through the levels is realized along the continu- 
ous paths of Fig.[2Ib) in the e, <f> space. The distribution 
function at each level is again constant though its devi- 
ation from equilibrium is the largest: For particles mov- 
ing down (up) the levels it coincides with the equilibrium 
1 - 2/„ = ±tanh(|A|/2T) at e = ±|A|. At this point, 
the dc current reaches its maximum. For a continuous 
path Eqs. ©, JTT} yield (3 + 7 = ±(<f> - tt)/2 + im and 



de/d4> = ±(hv x /2d)\S 12 \ 



(12) 



Transitions from one level to another at the crossing 
points <j) = irk are accompanied by sign changes of 
sin[2(/3 + 7)] in Eq. i|ll|) while the sign of de n /dcb is un- 
changed; in Eq. (|12|) it is plus for transitions upwards 
and minus otherwise. Since the variation in energy for 
each level is small compared to A we finally obtain from 
Eq. JjJJ replacing the sum with the integral 



^tanhf^ 
d \2T 

8e|A| T 2 L 



|A| 



|A| 



\S l2 (e)\ 2 ^de 



dc 



: tanh - — - 

Trh Vl-T 2 V2T 



(13) 



where L = In [cot(<$/2)]. The density of states is dn/de — 
n d/3/de = djithv x . 

Equation (|13|l accounts for the dc current from sub- 
gap levels (the corresponding superscript is omitted). 
Continuum states do also contribute. Indeed, for low 
transparency the Josephson current is 0, Q Ij = 
[ev x T 2 Y(a')/ird]smcbwhereY{a') = (2a'-7r)/ sin(2a'- 
7r) and < a' < it. Dc contributions of the type of 
Shapiro steps appear when the transparency T 2 Y(a') is 
ac modulated in resonance with loj. However, the mag- 
nitude of these steps is proportional to ev x T 2 /d and is 
thus much smaller than the dc current in Eq. i|13|) . 

Discussion.-lt is instructive to look at the energy bal- 
ance Eq. (JHJ. Since cb— ir = ujjt, the level sum of separate 
averages of the second term in the r.h.s. can be writ- 
ten as an integral over one continuous path e(<ft) (se e, 
Fig. [3(b) ) . For particles moving upwards, ^2 n de n /dt = 
(27T)" 1 f*+(de/dcp) dcj>, where e(0±) = ±|A|. Similarly to 
adiabatic pumps 0, Q this integral can be viewed as a 
circulation of a "vector potential" de n /dcb along a closed 
contour in the plane of complex A (circle of constant 
|A|). Each continuous path provides a nonzero circula- 
tion for an adiabatic process that takes place between 
two dissipative events when the path merges with con- 
tinuum. This term thus gives a contribution similar to 
that of multiple Andereev reflection processes in ballis- 
tic SNS junctions ^5|; in our case the barrier- induced 
minigaps are closed by the resonance gate voltage. The 



work by the gate voltage (first term) in Eq. (|SJ is an en- 
ergy counterflow; it reduces the pumped charge as com- 
pared to its ballistic value. Indeed, J2 n Q n (dV g /dt) = 
— (ithv x /ed) J2 n Q n - The average current is thus propor- 
tional to the Andreev probability, — (irhv x /ed) J2 n ( e + 
Qn) = —(ntiVx/d) J2 n l^ial 2 which agrees with Eq. JT3J. 

The energy varies from — |A| to +|A| during time Tq ~ 
I A\d/hvp luj which should be shorter than the inelastic 
relaxation time r e . Combined with the adiabatic condi- 
tion this gives the Josephson frequency window for the 
first resonance, \ A\d/hv F T e < 2eV/h < v F T 2 {l~T 2 )/d. 
Such slow r e _1 requires the use of low temperatures. 

In the case of a sinusoidal gate voltage Eq. QJ, res- 
onances also occur for rational gate voltage amplitudes 
specified by n and for luj = 2Q(1 + 2m). The dc current- 
voltage curve has peaks I m ~ / mal /(l + 2m) decay- 
ing with 77i but weakly depending on 77. The frequency 
2eV/h(l + 2m) should also belong to the window above. 
This puts an upper limit on observation of peaks of higher 
order in 777. Ideally, the largest peak I max is N times 
higher than the usual Shapiro step. 

To summarize, we found a new mechanism of nonequi- 
librium charge pumping. It is based on the spectral flow 
in a voltage-biased SINIS quantum junction when the 
minigaps in the energy spectrum are closed sequentially 
in resonance with the Josephson frequency. 
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